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ABSTRACT 

We wish to estimate magnetic field strengths of radio emitting galaxy clusters by minimising 
the non-thermal energy density contained in cosmic ray electrons (CRe), protons (CRp), and 
magnetic fields. The classical minimum energy estimate can be constructed independently of 
the origin of the radio synchrotron emitting CRe yielding thus an absolute minimum of the 
non-thermal energy density. Provided the observed synchrotron emission is generated by a 
CRe population originating from hadronic interactions of CRp with the ambient thermal gas 
of the intra-cluster medium, the parameter space of the classical scenario can be tightened by 
means of the hadronic minimum energy criterion. For both approaches, we derive the theoreti- 
cally expected tolerance regions for the inferred minimum energy densities. Application to the 
radio halo of the Coma cluster and the radio mini-halo of the Perseus cluster yields equipar- 
tition between cosmic rays and magnetic fields within the expected tolerance regions. In the 
hadronic scenario, the inferred central magnetic field strength ranges from 2.4 (Coma) to 
8.8 |J.G (Perseus), while the optimal CRp energy density is constrained to 2% ± 1% of the 
thermal energy density (Perseus). We discuss the possibility of a hadronic origin of the Coma 
radio halo while current observations favour such a scenario for the Perseus radio mini-halo. 
Combining future expected detections of radio synchrotron, hard X-ray inverse Compton, and 
hadronically induced y-ray emission should allow an estimate of volume averaged cluster 
magnetic fields and provide information about their dynamical state. 

Key words: magnetic fields, cosmic rays, radiation mechanisms: non-thermal, elementary 
particles, galaxies: cluster: individual: Coma (A 1656), Perseus (A 426) 



1 INTRODUCTION 

Clusters of galaxies harbour magnetised plasma. In particular, the 
detection of diffuse synchrotron radiation from radio halos or relics 
provides evidence for the existence of magnetic fields within the 
intra- cluster medium (ICM) (for a review, see ICarilli & Tavloil 
l2002h . Since the detection rate of radio halos in galaxy clusters 
seems to be of the or der of 30% for X-ray luminous clusters 
jGiovanninietalJll999t) . the presence of ma gnetic fields appears 
to be c ommon. Based on these observations. lEnfilin & Rottgerind 
llOO'i) developed a redshift dependent radio halo luminosity func- 
tion and predicted large numbers of radio halos to be detected with 
future radio telescopes. 

A different piece of evidence comes from Faraday rotation 
which arises owing to the birefringent property of magnetised 
plasma causing the plane of polarisation to rotate for a nonzero 
magnetic field component along the propagation direction of the 
photons tClarke et alj|200lh . However, the accessible finite win- 
dows given by the extent of the sources emitting polarised radi- 
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ation are a limitation of this method. The derived magnetic field 
strengths depend on the unknown magnetic field autocorrelation 
length which has to be deprojected from the observed two dimen- 
sional Faraday rotation measure maps using certain assumptions 
(see, however, EnBlin & Voet 2003; Vost & EnBlin 2003). A dif- 
ferent approach is given by the energy equipartition argument if a 
particular cluster exhibits diffuse radio synchrotron emission. The 
method assumes equal energy densities of cosmic ray electrons and 
magnetic fields in order to estimate volume averaged magnetic field 
strengths. 

The minimum energy criterion is a complementary method. It 
is based on the idea of minimising the non-thermal energy density 
contained in cosmic ray electrons (CRe), protons (CRp), and mag- 
netic fields by varying the magnetic field strength. As one boundary 
condition, the implied synchrotron emissivity is required to match 
the observed value. Additionally, a second boundary condition is 
required mathematically which couples CRp and CRe. For the clas- 
sical case, a constant scaling factor between CRp and CRe energy 
densities is assumed. However, if the physical connection between 
CRp and CRe is known or assumed, a physically better motivated 
criterion can be formulated. As such a case, we introduce the mini- 



2 C. Pfrommer and T. A. Enfilin 



mum energy criterion within tlie scenario of hadronically generated 
CRe. 

Classically, the equipartition/minimum energy formulae use 
a fixed interval in radio frequency in order to estimate the total 
energy density in cosmic ray electrons (CRe), a purely observa - 
tionally motivated procedure jBurbidgelll95a |Pacholczv]dll97Cl) . 
However, this approach has a drawback when comparing different 
field strengths between galaxy clusters because a given frequency 
interval corresponds to differen t CRe energ y intervals depending 
on the magnetic field strengths iBeckl2001l) . For this reason, vari- 
ants of the minimum energy criterion have been studied in order to 
place the magnetic field estimates on more physical grounds, based 
then on assump tions such as the fixed interval in CRe energy I Pohl 
ll993tlBeck et al. 1996; Brunetti et aL 1997') . The modified classi- 
cal minimum energy criterion does not specify a particular energy 
reservoir of the CRe. However, this apparent model-independence 
is bought dearly at the cost of the inferred magnetic field strength 
depending on unknown parameters like the lower energy cutoff of 
the CRe population or the unknown contribution of CRp to the non- 
thermal energy density. In the following, we use the term classical 
minimum energy criterion in its modified version, including e.g. a 
fixed interval in CRe energy. 

Natural candidates for acceleration mechanisms providing a 
highly-relativistic particle po pulation are strong structur e forma- 
tion and merger shocks fe.g.. lHarris et alj|l98CtlSarazin| [l999) or 
reacceleration processes of 'mildly' relativistic CRe {y^ ^ 100 - 
300) being injected over cosmological timescales into the ICM. 
Owing to their long lifetimes of a few times 10' yea rs, these 
mildl y relativistic CRe can accumulate within the ICM (see lSarazii] 
I2OO2I and references therein), until they experience continuous in- 
situ acceleration via resonant pitch angle scattering by turbulent 
Alfven waves as originally prop o sed bv Jaffe 1 1977.) and reconsid- 
ered b v lSchlickeiseret alHl987h.lBrunetti et alJj200lh . lo'hno et alj 
ilQQf). 'Oitti et al.' ^2002). and IKuo et alj J2003h . However, this 
reacceleration scenario also faces challenges as recent results im- 
ply: Isrunetti^talj |2004) show, that if the CRp-to-thermal energy 
density ratio were more than a few percent, Alfven waves would 
be damped efficiently such that the reacceleration mechanism of 
the electrons is inefficient. Because nearly all conceivable electron 
acceleration mechanisms produce a population of CRp which ac- 
cumulates within the clusters volu me, this re presents an efficient 
damping source of Alfven waves.' iKuo et aljT2004) presented an 
interesting line of argumentation to investigate the nature of radio 
halos by comparing the observed and statistically predicted popu- 
lation. This approach might allow to measure the life time of radio 
halos and thus help to conclude their physical origin with a future 
flux-limited, controlled, and homogeneous radio halo sample. 

In this work, we examine a minimum energy criterion within 
another specific model for the observed extended radio halos 
of ~ Mpc size: Hadronic interactions of CRp with the ambi- 
ent thermal gas produce secondary electrons, neutrinos, and 7- 
rays by inelastic collisions taking place throughout the clus- 
ter volume. These secondary CRe would generate radio halos 
through synchrotron emission (Dennison 198(7; 'Vestrand 1982; 
[ Blasi & Colafrancesco. .l999.; Dolag & En Blin 2000; MiniatLeXaL, 
boOlUPfrommer & Enfilirj2004 . This scenario is motivated by the 
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following argument: The radiative lifetime of a CRe population 
in the ICM, generated by direct shock acceleration, is of the or- 
der of 10** years for ~ \(f . This is relatively short compared 
to the required diffusion timescale needed to account for such ex- 
tended radio phenomena 1 Brun ettill2002h . On the other hand, the 
CRp are characterised by lifetimes of the order of the Hubble time, 
which is long enough to diffuse away from the production site and 
to be distributed throughout the cluster volume to which they are 
confined by magnetic fields dVolk et al. 1996; EnBlin eJ^aL 19971 
berezinskv et aljl997h . The magnetic field strength within this sce- 
nario is obtained by analogy with the classical minimum energy cri- 
terion while combining the CRp and CRe energy densities through 
their physically connecting process. Apart from relying on the par- 
ticular model, the inferred magnetic field strengths do not depend 
strongly on unknown parameters in this model. - 

The philosophy of these approaches is to provide a criterion 
for the energetically least expensive radio synchrotron emission 
model possible for a given physically motivated scenario. To our 
knowledge, there is no first principle enforcing this minimum state 
to be realised in Nature. However, our minimum energy estimates 
are interesting in two respects: First, these estimates allow scrutin- 
ising the hadronic model for extended radio synchrotron emission 
in clusters of galaxies. If it turns out that the required minimum 
non-thermal energy densities are too large compared to the thermal 
energy density, the hadronic scenario will become implausible to 
account for the extended diffuse radio emission. For the classical 
minimum energy estimate, such a comparison can yield constraints 
on the accessible parameter space spanned by lower energy cut- 
off of the CRe population or the contribution of CRp to the non- 
thermal energy density. Secondly, should the hadronic scenario be 
confirmed, the minimum energy estimates allow testing for the re- 
alisation of the minimum energy state for a given independent mea- 
surement of the magnetic field strength. 

This article is organised as follows: After introducing syn- 
chrotron radiation of CRe (Sect. l2.l'l . analytic formulae for hadron- 
ically induced emission processes are presented (Sect. 12. 2> . The 
classical and hadronic minimum energy criteria are then derived, 
the theoretically expected tolerance regions are given, and limit- 
ing cases are discussed (Sect.|^. In Sect. |31 we examine whether 
future observations of inverse Compton emission and hadronically 
induced y-ray emission can serve as tests for the verification of the 
minimum energy criterion. Magnetic and cosmic ray energy den- 
sities and their tolerance regions are inferred from application of 
the minimum energy arguments in the Coma and Perseus cluster 
for both scenarios (Sect. |5}. Our article concludes with formulae 
which provide recipes for estimating the magnetic field strength in 
typical observational situations (Sect.|Sl. Throughout this paper we 
use the present Hubble constant Hq = 70 h^Q km s"' Mpc"', where 
h^o indicates the scaling of Hq. 



2 THEORETICAL BACKGROUND 

This section presents our definitions and the theoretical background 
for this work. After introducing characteristics of the CRe popula- 
tion and the synchrotron emission formulae, we focus on specifi- 
cations of the CRp population. Finally, the section concludes with 

^ Likewise the minimum energy criterion within the reacceleration sce- 
nario of mifdfy refativistic CRe (ye - 100 - 300) can be obtained by min- 
imising the sum of magnetic, mifdfy refativistic CRe, and turbufent energy 
densities white affowing for constant synchrotron emission. 
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analytic formulae describing the hadronically induced y-ray and 
radio synchrotron emission processes. 

2.1 Cosmic ray electrons and synchrotron emission 

The differential number density distribution of a CRe population 
above a MeV is often represented by a power-law in energy E^, 



/e(r,£e)d£e dV = ncReW 



dV. 



(2.1) 
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where the tilde indicates that ncRc is not a real CRe number density 
although it exhibits the appropriate dimensions. The normalisation 
«CRe('') might be determined by assuming that the kinetic CRe en- 
ergy density ecRc('') is expressed in terms of the thermal energy 
density eth('')> 



ecRc(»') 



GeV 



VGeV/ 



Here we introduced the abbreviation 
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in order to account for cutoffs of the CRe population. If the CRe 
popul ation had tim e to lose energy by means of Coulomb interac- 
tions <GoulJl972h . the low energy part of the spectrum would be 
modified. This modification, which impacts on the CRe distribu- 
tion function f^(r,E^) and thus on A^^iae), can be approximately 
treated by imposing a time de pendent lower energy cutoff on the 
CRe population as described in Pfr ommer & En Blin 12004'). 

While the functional dependence of the CRe scaling parame- 
ter Xqrc (*■) is a priori unknown, its radial behaviour will be adjusted 
such that it obeys the minimum energy criterion. The thermal en- 
ergy density of the ICM e,), is given by 
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(2.5) 
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counts the number of particles per electron in the ICM using the 
primordial "'He mass fraction = 0.24. and denote the 
electron temperature and number density, respectively. 

The synchrotron emissivity jy at frequency v and per steradian 
of such a CRe population <2.I> . which is l ocated in an isotropic 
distribution of magnetic fields (Eq. (6.36) in lRvbicki & LightmanI 
Il979 1. is obtained after averaging over an isotropic distribution of 
electron pitch angles yielding 

n(ffv+l)/2 
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(2.9) 

(2.10) 



where r(a) denotes the Gamma-function ^Ab ramowitz & StegunI 
1 1 965ft and Oy - (o^e- 1)/2. Note that for later convenience, we intro- 
duce a (frequency dependent) characteristic magnetic field strength 
fic which implies a characteristic magnetic energy density Eg^. 
Line-of-sight integration of the radio emissivity jXr) yields the sur- 
face brightness of the radio emission 5,,(rx). 



2.2 Cosmic ray protons 

2.2.1 CRp population 

In contrast to the previously introduced CRe population and ow- 
ing to the higher rest mass of protons, we assume the differential 
number density distribution of a CRp population to be described 
by a power-law in momentum pp which for instance is motivated 
by shock acceleration studies: 

cdp^\ 



ff(r, pp) dpp dV = iicRfir) ^^-^ j 



GeV 



dV. 
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The normalisation ncRp('') can be determined in such a way that 
the kinetic CRp energy density £cRp('') is expressed in terms of the 
thermal energy density £th('') of the ICM, 



ecRp('') 



AeMp) 



XcRp(r)Ett,(r) = Ae^ ncRp(r) 
2(Q'p- 1) \ GeV 



- 2 3-a„ 



(2.12) 
(2.13) 



S{a,b) denotes the Beta-function I'Abramowitz & Stegun 1965). 

Aging imprints a modulation on the low energy part of the 
CRp spectrum by Coulomb losses in the plasma. This modification, 
which impacts on the CRp distribution function /p(r, pp) and thus 
on A£p(c!'p), can be treated appro ximately by imposing a lo wer mo- 
mentum cutoff as described in Pfrommer & EnBlii] i2004ft . On the 
other side, highly energetic CRp with energies beyond 2x 10' GeV 
are able to escape from the galaxy cluster assuming momen- 
tum dependent CRp difi"usion in a turbulent magnetic field with 
a Ko lmogorov-type spectrum on small scales ( Berezinskv et 4] 
Il997ft . The finite lifetime and size of particle accelerating shocks 
also give rise to high-energy breaks in the CRp spectrum. These 
low and high momentum cutoff's are always present in the CRp 
population. However, for CRp spectral indices between 2 < dp < 3 
these spectral breaks have negligible influence on the CRp energy 
density. If the breaks are neglected, the CRp energy density would 
diverge for ttp < 2 at the high-energy and for 3 < o-p at the low- 
energy part of the spectrum. In these cases, breaks have to be in- 
cluded by replacing Ae in Eq. )2.12t with 
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(2.14) 

(2.15) 

Beta-function 
are the lower 



where Bx{a,b) denotes the incomplete 
I Abramowitz & Stegu 3 Il96^ and p\ and 
and higher break momenta, respectively. 

2.2.2 Hadronically induced y-ray emission 

The CRp interact hadronically with the ambient thermal gas and 
produce pions, provided their momentum exceeds the kinematic 
threshold p,i,r = 0.78 GeV c"' for the reaction. The neutral pions 
decay into y-rays while the charged pions decay into secondary 
electrons (and neutrinos): 

n- p- + Vfj/v^ -> + vjv, -1- + 

n° ly. 

Only the CRp population above the kinematic threshold pd,, is vis- 
ible through its decay products in y-ray and synchrotron emission. 



4 C. Pfrommer and T. A. Enfilin 



An analytic formula describing the omnidirectional (i.e. inte- 
grated over 47t solid angle) differential y-ray source function re- 
sulting from /r^-decay of a power-law CRp population is given in 
pfro mmer & EnBlin (2004): 



qy{r, Ey) AEy iV CTpp c «n('') 2 

2Ey 



3 



IUjA c- 

GeV 

d, 



«CRpW 

GeV 

2Ey 



(2.16) 



dEy dV, 



-Xhc) denotes the target nu- 



where nt^(r) = 4^. n^ir) = n^{r)l{\ ^ 
cleon density in the ICM while assuming primordial element com- 
position with Xne = 0.24, which holds approximately. The formal- 
ism also includes the detailed physical processes at the threshold of 
pion production like the velocity distribution of CRp, momentum 
dependent inelastic CRp-p cross section, and kaon decay channels. 
The shape parameter 6y and the effective cross section cr^p depend 
on the spectral index of the y-ray spectrum ay according to 



0. 14 tt: 



' -I- 0.44 and 

4.4- 2.4 ff. 



(2.17) 

o-pp ^ 32-(0.96 + e''*-2"'">') mbam. (2.18) 

There is a detailed discussion in lPfrommer & EnfilinI i2004h how 
the y-ray spectral index ay relates to the spectral index of the par- 
ent CRp population ttp. In Dermer's model, which is motivated by 
accelerator experiments, the pion multiplicity is independent of en- 
ergy yielding the relation ay = a^ IPerme r 1986a b). 

Provided the CRp population has a power-law spectrum, the 
integrated y-ray source density Ay for pion decay induced y-rays 
can be obtained by integrating the y-ray source function qy{r, Ey) 
in Eq. <2.16> over an energy interval yielding 

I dEy qy(r, Ey) 

Jei 

Ayiap) A-'{ap) XcRp(r) njir) kT,(r), 



Ay(r, Ei,E2) 



(2.19) 
(2.20) 
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The y-ray number flux Ty is derived by means of volume integra- 
tion over the emission region and correct accounting for the growth 
of the area of the emission sphere on which the photons are dis- 
tributed: 



TyiEuE.): 



1 + z r 

AnD^ J 



dV Ay[r,(l+z)Eu{l+z)E2\. 
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Here D denotes the luminosity distance and the additional factors 
of I + z account for the cosmological redshift of the photons. 



2.2.3 Hadronically induced synchrotron emission 

Following for i nstanc e iDolag & EnlSlinl J2000l) and 
IPfrommer & Enfiiii] i2004h . the steady-state CRe spectrum is 
governed by injection of secondaries and cooling processes so that 
it can be described by the continuity equation 

[E,{r, E,)Ur, E,)) = q,(r, E,) . (2.24) 

For £'c(r, p) < this equation is solved by 



Ur,E,) ■■ 



1 



\E,{r,E,)\ 
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For the energy range of interest, the cooling of the radio emitting 
CRe is dominated by synchrotron and inverse Compton losses: 



4 (Tt c 7 
- E^{r, E^) = [Egir) + ecmb] El 

i mi C* 



(2.26) 



where is the Thomson cross section, SB{r) is the local magnetic 
field energy density, and ecmb = Scmb / energy density of 
the cosmic microwave background expressed by an equivalent field 
strength ficMB = 3.24(1 + z)"nG. Assuming that the parent CRp 
population is represented by a power-law )2.11> . the CRe popula- 
tion above a GeV is therefore described by a power-law spectrum 



Ur.E,) 
and hcK^ir) 



ncKAr) I E, Y' 
GeV \GeV/ ' 

eair) + Scmb 
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(2.27) 
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(2.29) 



a^-2 (Tt GeV 

where the effective CRp-p cross section cTpp is given by Eq. )2.181 . 

The hadronically induced synchrotron emissivity jy at fre- 
quency V and per steradian of such a CRe population <2.27t which 
is located in an isotropic distribution of magnetic fields within the 
halo volume is given by Eq. <2.7> . However, the normalisation 
hcReix) of the CRe population is given by Eq. <2.28> . The spectral 
index of the synchrotron emission is related to the CRp spectral 
index by a^ = (a^ - l)/2 = Q'p/2. 



3 MINIMUM ENERGY CRITERIA 

This section develops minimum energy criteria in order to esti- 
mate the magnetic field and studies the tolerance region of the ob- 
tained estimates. As described in Sect.Q we discuss two different 
approaches when requiring the non-thermal energy density of the 
source to be minimal for a particular (observed) synchrotron emis- 



The classical minimum energ 


y criterion used in radio astron 


Pohl 1993; Becketal. 1996; 


Brunetti et al. 1997) can be ap 



plied irrespective of the particular acceleration process of CRe, but 
unfortunately it relies on uncertain assumptions or parameters. Pro- 
vided the hadronic scenario of synchrotron emission applies, these 
dependencies can be softened. The resulting minimum energy ar- 
gument needs to be changed accordingly. 

The philosophy of these approaches is to provide an estimate 
of the energetically least expensive radio emission model possi- 
ble in each of these physically motivated scenarios. Thus, the ob- 
tained minimum energy estimates should not be taken literally in 
a sense that they are necessarily realised in Nature. However, the 
minimum energy estimates allow scrutinising the hadronic model 
for extended radio synchrotron emission in clusters of galaxies by 
comparing to the thermal energy density. For the classical mini- 
mum energy estimate, such a comparison can yield important con- 
straints on the accessible parameter space. 

The non-thermal energy density in the intra-cluster medium 
(ICM), which is the quantity to be minimised, is composed of the 
sum of the energy densities in magnetic fields, CRp and CRe: 



£nT - £s + EcRp + £cRc 



(3.1) 



The CRp population also includes higher mass nuclei in addition 
to protons. For convenience, we introduce canonical dimensionless 
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energy densities by means of scaling with the critical magnetic en- 
ergy density Ss^ = Bl/(Sn), where 6;. is defined in Eq. <2.9> : 

^ntC*") = — , -"(CRpir) = — , xcRc(r) = — , (3.2) 



, xth(r) = , and xcmb = 

^Bc ^Bc ^B^ 



(3.3) 



After presenting the conceptually simpler classical minimum en- 
ergy criterion, we will subsequently discuss the hadronic minimum 
energy criterion. 



The corresponding CRp, CRe and magnetic energy densities are 



given 


by 
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r nl/(2+i5) 


(3 


XcRe^ 
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-- C,„,,(r)[(l+fcp)(l-H^)C„a..(r)] 


-{l+S)/{2+S) 


XCRp,„ 






(3 



Note that these formulae deviate from lPacholczvld il970l) . since we 
use a fixed interval in CRe energy rather than in radio frequency 
(c.f. Sect, [p . Howe ver, they are equivalent to those obtained by 
lBrunettiet"a^^l997^ . 



3.1 Classical minimum energy criterion 

This section presents the classical minimum energy criterion from 
a physically motivated point of view, probing synchrotron emission 
from CRe characterised by a power-law distribution function with- 
out specifying their origin. This approach implies putting up with a 
dependence of the inferred magnetic field strength on unknown pa- 
rameters like the lower energy cutoff of the CRe population or the 
unknown contribution of CRp to the non-thermal energy density. 



3.1.1 Derivation 

Assuming a proportionality between the CRe and CRp energy den- 
sities, i.e. xcRp = ^pXcRe, the non-thermal energy composition equa- 
tion <3.H can be written as 



xnt - xs + (1 + fcn)XcRc 



(3.4) 



This assumption is reasonable if for instance the thermal electron 
population and the CRp were energised by the same shock wave 
assuming that there is a constant fraction of energy going into the 
CRp population by such an acceleration process, provided injection 
processes alone determine the energies. In order to proceed, we 
need an expression for jy (given by Eq. O) 

as a function of the 
dimensionless energy densities considered here: 

Jv(r) - xcRc{r)XB(r) 



S = 



- 1 



where 



(3.5) 
(3.6) 



2 4 

For consistency reasons, which will become clear in Sect. 13.21 we 
introduce <5 as a parametrisation of the spectral index. For typical 
radio (mini-)halos, S is a small quantity: The synchrotron spectral 
index Oy = 1 corresponds to S = 0. Solving Eq. <3.5> for the CRe 
energy density yields 



XcRc 



^£c Jv(r) jy 



where 



(3.7) 
(3.8) 



is a convenient auxiliary variable. Combining Eqs. 13. 4t and 13. 7> 
yields the non-thermal energy density solely as a function of the 
magnetic energy density 



XnT = Xfi + (1 + kp) Cdass(»') XB(r) 



(3.9) 



Requesting this energy density to be minimal for a given syn- 
chrotron emissivity yields the energetically least expensive radio 
emission model possible in this approach: 



9xnt 



9xb 



1 - (1 + k„)(l + ^)Cdass(/-)XB(r)-'-' = 0. 



(3.10) 



3.1.2 Localisation of classical minimum energy densities 

We wish to quantify how tight our statements about the inferred 
minimum energy densities are, i.e. to assign a tolerance region to 
the minimum energy estimates. This region would have the mean- 
ing of a quasi-optimal realisation of the particular energy densities. 
The curvature radius at the extremal value is one possible way of 
characterising the 'sharpness' of the minimum:^ 

\-l/2 



1 9~xnt 




^ xnt 9xg 


->"S„,;„ / 



(3.14) 



In order to avoid unphysical negative values for the lower tolerance 
level of xb^-^ we rather adopt the following logarithmic measure of 
the curvature: 



0"ln 



9^ In Xnt 




[9(lnxs)2 





-1/2 



(3.15) 



Considering the linear representation of Xb, this definition explic- 
itly implies tolerance levels which are given by exp(lnxB ± (ri„.vg). 
Applying this definition to Eq. <3.9> yields the theoretical tolerance 
level of the estimated minimum magnetic energy density. 



Xfi;;_'[Cd... (i + fcp) + x^;;__] 

Cdass(l+^p) (2 +.5)2 



2n 1/2 



(3.16) 



The tolerance level of the estimated minimum CRp energy density 
is given by 



0"ln 



9 In xcRp 



9 Inxg 



O" In.v 



(3.17) 



while the theoretical tolerance level of the estimated minimum CRe 
energy density can be obtained likewise. Applying this Gaussian 
error propagation to Eq. <3.7t . we obtain following general result: 



(Tin. 



CRp 



0"ln.vcRc = (1 + <5)0"ln.VB- 



(3.18) 



3.1.3 Equipartition condition 

In order to investigate under which conditions the classical mini- 
mum energy criterion implies exact equipartition, we examine the 



^ One could picture this approach by assuming Gaussian statistics for the 
distribution of non-thermal energy densities. The curvature radius at the 
minimum would then correspond to the width erg and thus yielding the 
68%-confidence level with respect to this extremal value. 
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special case (5 = of Eq. <3.1 1> . The resulting minimal magnetic 
and CRp energy densities read 



XB„.„„(r) = ^(1 +fcp)C,,assW, 
-*^CRpo.„,i„('') = 



(3.19) 
(3.20) 



This comparison shows that there exists exact equipartition be- 
tween the CRp and magnetic energy densitie s if S = and 
fcn 5S> 1! In our Galaxy ytp ^ 100 I Beck et ai]|l996h and ttv - 
ttGHz = 0.8 suggesting S = -0.1 iBeuer mann et alj 'l985). which 
implies that these equipartition conditions are well fulfilled. Com- 
paring this result with studies that are using a combination of syn- 
chrotron emission, the local CRe density, and diffuse continuum 
y-rays, [strong et al. 1 2000) interestingly imply the sam e magnetic 
field strength as inferred from equipartition arguments iBeck et all 
Il996l) . The corresponding minimal CRe energy density is given by 



-«"CRei 



(1+U ~^ 



(3.21) 



The classical minimum energy densities in CRe and magnetic fields 
are in exact equipartition if (5 = and kp <g: 1 ! In the limiting case 
6 = 0, the theoretical tolerance levels of the estimated classical 
minimum energy densities read 



■ CTln 



■VCRpo 



CTln 



1. 



(3.22) 



3.2 Hadronic minimum energy criterion 



This section deals with the hadronic minimum energy criterion and 
thus probes hadronically induced synchrotron emission. If a sig- 
nificant part of the CRe population is known to be generated by 
hadronic interactions, the minimum energy estimate allows testing 
for the realisation of the minimum energy state for a given inde- 
pendent measurement of the magnetic field strength. Thus, the pro- 
posed minimum energy criterion would provide information about 
the dynamical state of the magnetic field. 

3.2.1 Derivation 

To derive the hadronic minimum energy criterion, we need an ex- 
pression for the CRp energy density as a function of the mag- 
netic energy density which is obtained by combining Eqs. <2.7> and 



Ch-idr{r) [xs(r) -I- xcmb] xgir) 
jv(r) jy(r) 



xcRpir) 



8 = 



ttv - 1 



a. 



-2 



(3.23) 
(3.24) 

The parameter Chadr('') has the meaning of a hadronic synchrotron 
emissivity per target nucleon density and per frequency. Its value 
decreases for any existing cutoff in the parent CRp population as 
described by Eq. J2.14t . For convenience we introduce the parame- 
ter (5 which ranges within [-0.1,0.2] for conceivable CRp spectral 
indices o-p = [1.6,2.8]. 

Assuming the hadronic electron source to be dominant, we 
can neglect the primary CRe population. In any case, the energy 
density of hadronically generated CRe is negligible compared to 
the energy density of the parent CRp population: Above energies 
of ~ GeV the differential hadronic kp = ncRp/"CRc has typically 



0.35 
0.30 
^ 0.25 



numerical solution, C|„dr 
asymptotic solution, Chadr = 1 
numerical solution, Chad, = 10"'' 
asymptotic solution, C|,adr = lO""* 




-0.10 -0.05 -0.00 0.05 

expansion parameter 6 = 



0.10 0.15 

(ffp - 2)/4 



Figure 1. Comparison between the numerical solution and the second order 
asymptotic expansion for Chadr = 1 (cool core cluster, formerly referred to 
as cooling flow cluster) and Chadr = lO"* (cluster without a cool core). The 
perfect match indicates a fast convergence of the asymptotic solution for the 
range of S being considered. 



values ranging between kp ~ 100 (Perseus) and kp ~ 300 (Coma), 
where we inserted the typical values of the central density and mag- 
netic field strength. The tilde in kp indicates the slightly modified 
definition compared to the previous classical case. Requiring min- 
imum non-thermal energy density within the hadronic framework 
for a given synchrotron emissivity yields 



9jcnt 



9xs 



= 1 + 



QjCCRp 



= 0. 



(3.26) 



By using Eq. <3.23> we obtain the following implicit equation for 
the minimum magnetic energy density: 

[d + 5) xcMB + SxB„^Jr)] XB^^Jr)-'-' = C^'^Jr). (3.27) 

The definition of Chadi('') <3.24> reveals an implicit dependence on 
the parametrised spectral index S. However, the right hand side of 
the minimum energy criterion <3.27> is uniquely determined for 
a given spectral index and an observed synchrotron emissivity at 
a particular frequency. Thus, the minimum energy density of the 
magnetic field giving rise to an observed synchrotron emission in 
the hadronic model can either be obtained by solving Eq. <3.27> 
numerically or applying the asymptotic expansion which will be 
developed in Sect. l3.2.2l 



3.2.2 Asymptotic expansion for 6 

The asymptotic expansion of the magnetic field energy density as a 
function of the small parameter S follows from the minimum energy 
criterion <3.27>: 



XBiir) 
XB.(r) 



XB„{r) + Sxb, (r) + S-XB,{r) + 0(8^ 



XBo I Xbi 

— \l-\nxB„ + 



XCMB 



(r). 



2 



-{lnxB,f -^{2 + \nxB,) + ^ 



(3.28) 
(3.29) 

(3.30) 
(3.31) 



A comparison between the numerical solution and the asymptotic 
solution for two different values of the parameter Chadr is shown 
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in Fig. Q The particular values of the parameter Chadr will be mo- 
tivated in Sect. 15.21 This figure should serve as mathematical il- 
lustration of the convergence behaviour of the asymptotic expan- 
sion while keeping the parameter Chadr fixed. The second order 
asymptotic solution perfectly agrees with the exact solution for 
the range 6 = [-0.15,0.2] of conceivable CRp spectral indices 
ttp = [1.6,2.8]. 

3.2.3 Localisation ofhadronic minimum energy densities 

Considering the accuracy of the estimated minimum energy den- 
sities in the hadronic scenario, we can also apply the logarithmic 
measure of the theoretical tolerance level as defined in Eq. <3.15> : 

0"ln XB (r) = [-4^1 + Chadr {xb„„„ + ^CMB )] (Chadr ^B^,„ 

X {Chadr-vcMB + ^1 [d + '5)'-^B„n + (2 + Sf xcmb])'"^ (r). (3.32) 

The corresponding tolerance level for the minimum CRp energy 
density is obtained by applying Gaussian error propagation <3.17> : 



, , . ^CMB 

0"ln.vcRp('') = 7T 

\Xb . (r) + XcMB 



+ S\ cri„,.g(j-). 



(3.33) 



The consequences of these rather unwieldy formulae will become 
intuitively clear in the next section where characteristic limiting 
cases are investigated. 



3.2.4 Special cases 

In order to gain insight into the hadronic minimum energy crite- 
rion, we investigate special cases of Eq. <3.27t . namely (5 — > 0, 
fig » EcMB, and Bb ■« Bcmb, while simultaneously considering the 
resulting tolerance regions of the previous Sect. l3.2.3l 

(i) (5 — > 0: The limit of (5 — > corresponds to a hard spectral CRp 
population described by a spectral index of ffp = 2.0 and thus o-,, = 
1 . There the dimensionless magnetic field energy density reads 



„('■) = VChadr('')-»''CMB ' 



jv(r) ECMB 



vnyiir) 



1/2 



(3.34) 



We can formulate the corresponding dimensionless CRp energy 
density resulting from this minimising argument: 



^CRp, 



„(»■) = Chadr(»') + VChadr(»')-«CMB- 



(3.35) 



It is interesting to note that there are two regimes for xcRp„^^^^, 
namely 



^CRp, 



„('•) 



Chadr('')> 



for Chadr(»-)» 0.01 



■^'Bo,min('')' for C'hadr('') 0.01 V^^^ . 



(3.36) 



Only in the limit of a small parameter Chadr('') and (5 = 0, there 
is an exact equipartition between the hadronic minimum energy 
densities in CRp and magnetic fields! 

The tolerance regions of the estimated minimum energy densities 
are also more intuitive to understand in the limit ^ = compared to 
the general case laid down in Eq. <3.32t : 



0"ln.VB„(»-) 



1 + ■ 



xb„ 



2j^cmb 



1/2 



1, 



2jt"CMB 



1/2 



implying two regimes: 

forXB„_^._, « XCMB, 
for Xb„^.^ » XcMB • 



(3.37) 



In the previously studied classical case <3.22t . the tolerance region 



remains constant for all conceivable magnetic energy densities. In 
contrast to this, the hadronic scenario shows a increasing toler- 
ance region cri„,.^^ for strong magnetic field strengths (fio.min ^ 
3.24 (1 -I- z)' nG) which is explained by the following argument: In 
the limit of strong magnetic fields, the hadronically induced syn- 
chrotron emission with a^, = 1 does not depend any more on the 
magnetic field strength but only on the CRp energy density. This is, 
because the inverse Compton cooling is negligible in this regime, 
implying that observed synchrotron emission is insensitive to any 
variation of the magnetic field strength since all injected CRe en- 
ergy results in synchrotron emission. Therefore magnetic field es- 
timates inferred from minimum energy arguments are rather uncer- 
tain in the limit of strong magnetic field strengths. 

The tolerance levels of the corresponding CRp energy density, 
derived from Eq. <3.33L shows two limiting regimes: 



0"lii 



O-lr 



TCRpo 



JtCRpo 



(r) 



1 + 



XCMB 



1 + 



^S„,,„i„('') 



1/2 



XcMB 



1/2 



2xcMB 
for -"^So .min ^ ^CMB 
for Xs„^,__ » XcMB 



implying 



(3.38) 



The tolerance region of the inferred minimum CRp energy den- 
sity decreases in the regime of strong magnetic fields which can 
be understood by the same token as above, i.e. synchrotron losses 
dominate over the inverse Compton cooling in this limit. Thus, the 
observed radio emission reflects accurately the CRp energy density 
in the strong magnetic field limit. 

(ii) Eg » Ecmb'- In the limit of » ^cmb^ the magnetic field 
energy density is an even stronger function of the synchrotron emis- 
sivity and the number density of the ambient gas: 



-^B„,„(r) 



[<5Chadr(r)]'/<'^'' 
^Chadr('-)'«'"^V 



jy(r) 



vnf^(r) 



1/(1+<S) 



(3.39a) 
(3.39b) 



where we assumed \S\ <c 1 in the second step. In the limit of strong 
magnetic fields and small 6, exact equipartition of the magnetic and 
CRp hadronic minimum energy density does not occur, because 

XCRp^Jr) = <5-'/<'^'"Chadr('-)"<'+''> - Chadr('-)""^'"(1 In5).(3.40) 

(iii) £cMB ^ eb'. In the opposite limit, we obtain the following 
minimum energy criterion for the magnetic field energy density: 



[(1 -(- 5) XcMB Chadr('')] 



1/(2+5) 



{'4) 



[^CMBChadr('')] 



1/(2+4) 



7v('')SCMB 



1/(2+5) 



(3.41a) 
(3.41b) 

(3.41c) 



vnt^(r) 

where we again assumed ^ <K 1 in the second step. In contrast to 
the previous limit, there is exact equipartition between the magnetic 
and CRp hadronic minimum energy density to zeroth order in d, 

(l+5)-<'^^>«^^*>[XcMBChadr(r)]'«^^^' 



|l ~ 2) [-''CMBChadr('')] 



1/(2+5) 



(3.42a) 
(3.42b) 



4 FUTURE TESTING 

This section will discuss possibilities of measuring the magnetic 
field strength, averaged over the cluster volume, in order to test for 
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the realisation of the energetically least expensive state given by the 
minimum energy criterion. 

4.1 Inverse Compton emission 

The CRe population seen in the radio band via synchrotron emis- 
sion should also scatter photons of the cosmic microwave back- 
ground (CMB), the local radiation field of elliptical galaxies, and 
the thermal X-ray e mission of t he IC M to different energy bands 
jpelten & Morrisor] |l966: Reej ligGTi) . Combining measurements 
of inverse Compton (IC) and synchrotron emission eliminates the 
uncertainty in number density of the CRe population provided the 
inevitable extrapolation of the CRe power-law distribution for cer- 
tain observed wavebands is justified. This enables the determina- 
tion of the magnetic field strength B for an IC detection and a lower 
limit on B for a given non-detection of the IC emission. 

The source function qic owing to IC scattering of CMB pho- 
tons ofi^ an isotropic powe r law distribution of CRe ( Eg. <2.1» is 
(derived from Eq. (7.31) in lRvbicki & Lightmanlig?^ in the case 
of Thomson scattering), 



q\c{r,Ey) 



q{r)f\c{ae) 



2 \ 1-ffc 



GeV; \^rcMB 
^ + 4a, + 11) 



-(ay+l) 



(tte -I- 3)2 (a, + 5){a, + l) 
+ 5\ ^1 a, + 5\ 



X r 



and q(r) = 



37to-T«cRi;('') (kTcua) 



(4.1) 



(4.2) 



(4.3) 



where Oy = (a^ — l)/2 denotes the spectra l index, f(fl) the Riemann 
^-function iAbramowitz & Stegunlll963) . and WcrcC'') is given by 
Eq. <2.2> . After integrating over the considered energy interval and 
the IC emitting volume in the cluster, the particle flux TyiEi , £2) is 
obtained (c.f. Eq. <2.23» . 

EnBlin & Biermann 1 1998) compiled non-detection limits of 
IC emission of different photon fields in various wavebands from 
the Coma cluster and obtained the tightest limits on B from 
the CMB photon field. The same CRe population emitting radio 
synchrotron radiation scatters CMB photons into the hard X-ray 
band. Non-detection of this IC emission by the OSSE experiment 
jRephaeU et"anil994h yields a lower limit on the central magnetic 
field strength of ficoma(O) > 0.2 \iG fg°*\ where /g is the filling 
factor of the magnetic field in the volume occupied by CRe. Pro- 
vided the CRe power-law distribution ca n be extrapo lated to lower 
energies, the limit given by the EUV flux iLieu et 31X1 996^ predicts 
a magnetic field strength stronger than 5coma(0) > 1.2 [lG /g" *''. 

The reported high ener gy X-ray excess of the Com a cluster 
by the Beppo-Sax satellite Ipusco-Pemiano et alJll99^ initiated 
other theoretical explanations about the origin of such an excess 
than IC up-scattering of CMB photons by relativistic electrons. 
One possibility implies the existence of a bremsstrahlung emit- 
ting supra-thermal electron population between 10 and 100 keV 
which would a lso produc e a unique Sunyaev-Zel'dovich signature 
iEnfilin et alJl999;: .EnBlin & Kais ea2000^.Blasi et al..2000^.Blasi> 
bOOdlLiang et all2002t la)lafrance sc'o et alj2003[l However, such 



a population is questioned on theoretical reasons iPetrosianl200lh . 
and even the high X-ray exc ess of Coma itself is un der debate 
jRossetti & Molendi 2004; Fu sco-Femiano et aljEo04l) . The data 
analysis of the RXTE observation of A 2256 yielded also evidence 
for a second spectral component (Rechaeli & Gruber 2003). On the 
basis of statistics alone, the detected emission is inconclusive as to 



whether it originates from a thermal multi-temperature fluid or an 
isothermal gas in combination with a non-thermal IC power-law 
emission. Future measurements with the IBIS instrument on-board 
the INTEGRAL satellite should provide even tighter upper limits 
respectively detections of the IC X-ray flux of a particular cluster 
and should therefore allow even tighter lower limits on the mag- 
netic field strength. 



4.2 7-ray emission 

This subsection outlines the method for estimating upper limits on 
the magnetic field strength using hadronic CRp interactions. The 
method is based on the idea of combining hadronically induced 7- 
ray and synchrotron emission to eliminate the uncertainty in num- 
ber density of the CRp population. For this purpose, one neces- 
sarily needs to resolve the detailed broad spectral signature of 7- 
rays resulting from the Tr^-decay (rt'^-bump centred on m^c-/2 ^ 
67.5 MeV) as laid down in Eq. <2.20t . This is to exclude other pos- 
sible processes contributing to diffuse extended 7-ray emission like 
IC radiation or dark matter annihilation. Because of possible other 
additional contributions to the diffuse synchrotron emission from 
CRe populations, e.g. primarily accelerated electrons, we are only 
able to provide an upper limit on the magnetic field strength. 

The proposed algorithm allows for different spatial resolutions 
of the y-ray and synchrotron emission. The application, we have in 
mind, is the determination of intracluster magnetic fields. In this 
case, y-ray observations of the /r^-decay induced 7-ray emission 
signature are only able to provide integrated 7-ray fluxes of the 
entire cluster due to their comparably large point spread function. 
7-ray fluxes depend on the thermal electron density and temper- 
ature profiles which have to be derived from X-ray observations. 
However, if we assumed a comparable resolution in y-ray and syn- 
chrotron emission the dependences on the thermal electron popula- 
tion could be eliminated (c.f. Sect. l2.2t . The algorithm consists of 
the following two steps: 

(i) Choosing a constant scaling parameter XcRp for the CRp pop- 
ulation and performing the volume integral of the energy integrated 
7-ray source density Ay <2.20> yields the y-ray flux according to 
Eq. <2.231 . The CRp parameter XcRp is obtained by comparing the 
observed to the theoretically expected y-ray flux. 

(ii) Inserting XcRp into the synchrotron emissivity jV <2.7t en- 
ables us to solve for the magnetic field strength as function of an- 
gle on the sky when comparing to radio surface brightness obser- 
vations. 

Once a detailed angular distribution of -decay induced 7- 
ray emission from a particular astrophysical object is available this 
algorithm may be implemented for the average of pixels contained 
within a certain solid angle. In this case the spatial distribution of 
CRp may even be deprojected. 

Is there a chance to apply this method to galaxy cluster mag- 
netic fields with future y-ray instruments? Because of the necessity 
of resolving the broad spectral signature of y-rays resulting from 
the ;r''-decay centred on ~ 67.5 MeV, the imaging atmospheric 
Cerenkov technique with a lower energy cutoff above 10 GeV is 
not applicable. Contrarily, the LAT instrument on-board GLAST 
scheduled to be launched in 2007 has an angular resolution better 
than 3.5° at 100 MeV while covering an energy range from 20 MeV 
up to 300 GeV with an energy resolution better than 10%. Assum- 
ing a photon spectral index of = 2 for the y-ray background, 
the point-source sensitivity at high galactic latitude in a one year 
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Table 1. Individual parameters describing the extended diffuse radio emis- 
sion in the Coma and Perseus galaxy cluster according to Eq. 15. 11 . The 
maximal radius to which these profiles are applicable is denoted by r„,ax. 
The radio data at 1.4 GHz are taken from Deiss et al. (1997) (Coma) 
and Pedlar et al. 1 1 990 ) (Perseus) w hile the profile of the Coma cluster at 
326 MHz is taken f romlGov orii et aU 1,2001.) which is based on radio obser- 
vations bv lVenturi et alJ Il990t) . 



(5.2) 



So 

Cluster [Jy arcmin"^] 


[ft^Q kpc] 


' max 

[/i-(J Mpc] 


/3 


1.4 GHz observations: 








A 1656 (Coma) 1.1x10"' 


450 


1.0 


0.78 


A426 (Perseus) 2.3 X 10"' 


30 


0.1 


0.55 


326 MHz observation: 

A1656 (Coma) 4.7 X lO"' 


850 


0.7 


1.07 



all-sky survey is better than 6 x 10"' cm"- s"' for energies inte- 
grated above 100 MeV. Specifically, assuming a CRp spectral in- 
dex Qfp = 2.3 and a flat Xcr^ for simplicity, such a one year all-sky 
survey is able to constrain XcRp < 0.01 (Perseus) and XcRp < 0.04 
(Coma). Taking additionally into account the y-ray flux between 
20 MeV and 100 MeV as well as a longer survey time will improve 
the sensitivities and yield even tighter limits on XcRp- Comparing 
these limits with energetically favoured values of ^cRp,,,;,, which 
are obtained by applying hadronic minimum energy arguments to 
a given radio synchrotron emission (c.f. Fig.|5J yields compara- 
ble values in the case of Perseus while the situation in Coma is 
less optimistic. However, a definitive answer to the applicability of 
this method can not be given on the basis of minimum energy ar- 
guments because such a minimum energy state is not necessarily 
realised in Nature. 



5 APPLICATIONS 

In this section we apply the classical and hadronic minimum en- 
ergy criterion to the radio (mini-)halos of the Coma and Perseus 
galaxy clusters. For simplicity, the CRp and CRe spectral indices 
are assumed to be independent of position and therefore constant 
over the cluster volume. If a radial spectral steepening as reported 
byloiovannini et al. ( 1993) in the case of the Coma radio halo will 
be confirmed by future radio observations evincing a better signal- 
to-noise ratio, the CRp and CRe spectral index distributions would 
have to embody this additional degree of freedom. We discuss in 
Sect. l5.3l that a moderate steepening would not significantly modify 
the hadronic minimum energy condition while a strong steepening 
would challenge the hadronic scenario. 



The individual parameters for the Coma radio halo and Perseus 
radio mini-halo are shown in Tabled Both profiles describe the 
extended diffuse emission where all point sources have been sub- 
tracted. Particularly, the extremely bright flat-spectrum core owing 
to relativistic outflows of the radio galaxy NGC 1275 in the centre 
of the Perseus cluster has been excluded from the fit. The electron 
density profiles are inferred from X -ray observations by Briel et ad] 
1 1992) (Coma) and lChurazov et aljj200.3t) (Perseus). 

In the following, we assume that the energy distribution of 
the CRe population above a MeV is represented by a power-law 
in energy with a lower cutoff. As a word of caution, such an 
assumption might be a strong simplification in the case of turbu- 
lent acceleration models yielding more complex e nergy distribu- 
tions which are considerabl y flatter at low er energies iBrunetti et alJ 
I2OOI; Petrosian 2001; Ohn o et aj2002h . We assume a CRe spec- 
tral index of o-e = 3.3 which translates into a synchrotron spectral 
index of cc,, = 1.15. This is consistent with radio data of Perseus and 
Coma; particularly when considering the spectral cutoff between 
1 and 10 GHz owing to th e Sunyaev-Zel'dovich flux decrement 
(Deis s et alll997l;lEnfilirl2 002; Pfrommer & EnBhn 2004). 

Applying the classical minimum energy criterion to the diffuse 
synchrotron emission of the Coma cluster yields a central magnetic 
field strength of 6coma('' = 0) = 1.1^0 4 l-'Gr. In the case of the 
Perseus radio mini-halo we obtain fipcrscus('' = 0) = V.2^2 g '^^^ 
indicated tolerance levels derive from the logarithmic definition of 
the theoretical accuracies of the minimum in Eq. <3.15t . However, 
these values are highly dependent on the lower energy cutoff of the 
CRe population, Ei , and the CRp proportionality parameter kp. Fol- 
lowing the philosophy of our paper we adopt a physically motivated 
lower Coulomb cutoff of the CRe distribution of Ei = 0. 1 GeV cor- 
responding to a relativistic 7 factor of ^ 200 (as suggested by 
laarazi 3l99i). We also adopt a conservative choice of the propor- 
tionality constant between the CRe and CRp energy densities of 
fcp = 1 . An increase of kp would directly increase the magnetic field 
strength by approximately the square root of this factor. 

For a cluster-wide comparison of energy densities of magnetic 
fields, CRe, and CRp, it is convenient to introduce a scaling with 
the thermal energy density by means of 



^sir) = — — ■ = — — ■, and 
x±{r) Etbir) 

/ \ _ -''CRp,CRe('') eCRp,CRi;('') 
ACRp.CRc('') = 7- - 7- . 

Xthir) Sthir) 



(5.3) 



(5.4) 



5.1 Classical minimum energy criterion 

In a first step we have to deproject the radio surface brightness and 
electron density profiles; In analogy to X-ray observations we as- 
sume the azimuthally averaged radio profile to be described by a 
j8-model. 



Sy(rj_) = So 



1 + 1^ 



-3/3+1/2 



(5.1) 



D eprojecting this profile yiel ds the radio emissivity (c.f. Appendix 
of IPfrommer & EnBliij20oi) 



Profiles of the CRe-to-thermal energy density XcRc^^^(r) and 
magnetic-to-thermal energy density Xs^^-^ir) are shown in Fig. |2| 
for the Coma and Perseus cluster and our different scenarios. The 
upper panels show the scaled energy densities in the acceleration 
model of CRe, obtained by the classical minimum energy crite- 
rion. While the optimal magnetic energy density is roughly a fac- 
tor of two larger than the CRe energy density, they both can be in 
equipartition for the quasi-optimal case of their distribution of en- 
ergy densities, as indicated by the dark shaded regions. In order to 
explain the observed synchrotron emission in the CRe acceleration 
scenario the CRe and magnetic energy densities are only required 
to be below one percent of the thermal energy density. 
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Coma cluster: classical minimum energy criterion 
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Perseus cluster: classical minimum energy criterion 
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Coma cluster: hadronic minimum energy criterion 



Perseus cluster: hadronic minimum energy criterion 
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Figure 2. Profiles of the CRe-to-thermal energy density Xcrc,^|„ (r) (solid) and magnetic-to-thermal energy density Xg^^.^ (r) (dotted) as a function of deprojected 
radius are shown. The different energy densities are obtained by means of the classical minimum energy ciiterion (upper panels) and the hadronic minimum 
energy criterion (lower panels). In the latter scenario, profiles of the scaled CRp energy density are shown instead of CRe profiles. The left hand side shows 
profiles of the Coma cluster while the right hand side represents profiles of the Perseus cluster The light shaded areas represent the logarithmic tolerance regions 
of XB„,i„('') and XcRe„,i„(r), respectively, while the dark shaded regions indicate the overlap and thus the possible equipartition regions in the quasi-optimal 



5.2 Hadronic minimum energy criterion 

Assuming a CRp spectral index of o-p = 2.3 when applying 
the hadronic minimum energy criterion to the diffuse synchrotron 
emission of the Coma cluster yields a central magnetic field 
strength of Scoma('' = 0) = 2.4^| p \iG. In the case of the Perseus 
cluster we obtain Bperseus('' = 0) = S.S;^^^^^ \iG. Both inferred pro- 
files of the magnetic field are relatively flat: While the magnetic 
field strength in the outer part of the radio mini-halo in Perseus 
{r ^ 100 /17Q kpc) declines to a value of 55% of its central value, 
the magnetic field in the outer region of the radio halo in Coma 
(r 1 h^^ IVIpc) only decreases to 72% of its central value. 



As discussed in Sect. l3.2.H the hadronic scenario shows an 
increasing tolerance region for strong magnetic field strengths 
which are expected to be present in the case of cool core clus- 
ters, such a s Perseus (Tavlor & Perlev 1993; Carilli & Taylor 20oS 
IVogt & En filin 2003j). In this limit, synchrotron losses dominate 
over inverse Compton cooling. This almost cancels the dependence 
of the synchrotron emissivity on the magnetic energy density. The 
lower panels of Fig.|2|show the scaled energy densities Xs^^Jr) and 



^CRp„,i„('') as inferred from the hadronic scenario. In both clusters 
the optimal CRp energy density is larger than the magnetic energy 
density within this model. However, both energy densities can be 
again in equipartition for the quasi-optimal case of their distribu- 
tion, as indicated by the dark shaded regions. 

Owing to the inverse dependence of cri„,.g and cr^^ ^^^^ on 
the magnetic energy density in the limit of strong magnetic fields 
(c.f. Eqs. <3.37t and <3.38» . a large tolerance region of Xb^-^ imme- 
diately implies a well defined localisation of -X'cRp^^.^ . In the Perseus 
cluster this results in a confinement for the CRp energy density of 
2% ± 1 % of the thermal energy density. On the other hand, in the 
Coma cluster Xg^„;„ {r) and ^cRp,,,,, ('') are required to increase by less 
than one order of magnitude from the centre to the outer parts of the 
cluster in order to account for the observed radio halo. This increase 
might be partly due to azimuthally averagin g the as pheric elec- 
trori density distribution of the Coma cluster jPfromm er & Enfil^ 

High values for the radio emissivity per target density and 
frequency Chadr of order unity seem to reflect conditions in cool 
core clusters (formerly referred to as cooling flow cluster) whereas 
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smaller values seem to represent conditions in clusters without cool 
cores as shown in the following: 

Chadr = -; — ~~r~ ~ Cclustci ( , . ) ( ~^ | ( I ' ^^'^^ 

AE,y„A^,a \ 1.4 GHz/ \;je.o/ \jv,ol 

where Ccoma = 9.4 x 10"" and Cpe,scus = 1-5 x 10"'. 

5.3 Possibility of a hadronic scenario in Perseus and Coma 

5.3.1 Perseus radio mini-halo 

The azimuthally averaged radio surface brightness profile of the 
Perseus mini-halo matches the expected emission by the hadronic 
scenario well on all radii (.Pfrommer & EnBlin.20Q4i) while requir- 
ing almost flat profiles for CRp and magnetic energy densities rela- 
tive to the thermal energy density, XcRp and Xb, respectively. More- 
over, the small amount of required energy density in cosmic ray 
protons £cRp (~ 2% relative to the thermal energy density) supports 
the hypothesis of a hadronic origin of the Perseus radio mini-halo 
not only because the hadronic minimum energy criterion predicts a 
close confinement of £cRp (see Sect. 15. 2t but also because cosmo- 
logical simulations carried out by Miniati et al. (2001) easily pre- 
dict a CRp population at the clusters centre of this order of magni- 
tude. 

5.3.2 Coma radio halo 

The energetically favoured radial profile for the magnetic field 
strength in the Coma cluster is almost flat as predicted by the 
hadronic minimum energy criterion (see Sect. 15. 2> . Provided these 
results would be realised in Nature, this apparently contradicts pro- 
files of the magnetic field strength as inferred from numerical sim- 
ulations which seem to follow the elec tron density ndr) accord ing 
to B{r) oc n^ir)"" with org 6 [0.5,0.9] fo olag et aljl 999'.'200n. It 
would also contradict theoretical considerations assuming the mag- 
netic field to b e frozen into the flow and isotropised, i.e. = 2/3 
jTribbldll993h . Applying the flux freezing conditions to the elec- 
tron density profile of Coma jBriel etalll 992) yields an expected 
decline of the magnetic field strength from its central value to the 
magnitude at 1 hj^^ Mpc by a factor of ~ 6.7. 

However, there are other numerical, physical, and observa- 
tional arguments indicating large uncertainties in the origin, am- 
plification mechanism, and specific profile of the magnetic field 
strength, thus leaving the hadronic scenario as a viable explana- 
tion of the Coma radio halo: In contrast to the cited numerical sim- 
ulations, there are other cosmological simulations <Miniatill200ll: 
iMiniati et al]l200lh which are able to produce giant radio halos 
in the hadronic scenario and therefore reasonably flat profiles of 
the magnetic field strength. From the physical point of view, there 
could be stronger shear flows or a larger number of weaker shocks 
in the outer parts of clusters which are unresolved or not accounted 
for in current simulations. This would imply stronger additional 
amplification of the magnetic field strength in the outer parts of 
clusters yielding a flatter profile of the magnetic field strength. Ob- 
servationally, there are still uncertainties in the radio surface pro- 
files which are increased by azimuthally averaging the diffuse syn- 
chrotron emission in the presence of non-centrally symmetric emis- 
sion components such as the so-called radio-bridge in Coma around 
NGC 4839. 

In order to account for the radio halo of Coma in the hadronic 
scenario, the product of XcRp and Xg needs to increase by nearly 
two orders of magnitude towards the outskirts of the halo (c.f. 
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Figure 3. Parameter study on the ability of hadronically originating CRe 
to account for the radio halo of Coma. Assuming the profile of the mag- 
netic field to scale with the square root of the electron density yields a flat 
magnetic-to-thermal energy density ratio Xb- Shown are contour lines from 
the bottom to the top of max(XcRp) = (1,0.3,0.1,0.03,0.01) for the range 
r < 1 hlj^ Mpc in parameter space spanned by o-p and Bq. Conservative 
choices for CRp spectral breaks have been assumed. The lower part rep- 
resents the region in parameter space, where the hadronic scenario faces 
serious challenges for explaining the observed radio halo of Coma. 

lower left panel of Fig. |2}. Leaving aside the minimum energy 
criterion, this increase can be split arbitrarily among the magnetic 
and CRp energy density ratios. For instance a constant magnetic- 
to-thermal energy density ratio Xb, corresponding to cfg = 0.5 
in an isothermal cluster, is still consistent within the theoretically 
expected tolerance regions, i.e. within a quasi-optimal realisation. 
However, the CRp-to-thermal energy density ratio XcRp would have 
to compensate for this by increasing nearly two orders of magni- 
tude towards the outskirts of the halo.'' This choice of the magnetic 
field morphology (o-g = 0.5) has been adopted in Fig. which 
represents a parameter study on the ability of hadronically origi- 
nating CRe to generate the radio halo of Coma. Contour lines of 
max(XcRp) = (1,0.3,0.1,0.03,0.01) for the range r < Ih'^ Mpc 
are shown in parameter space spanned by o-p and Bq. The gradient 
of the maximum of XcRp points downwards in Fig.|3]and thus leaves 
the upper region of paraineter space where the hadronic scenario is 
energetically able to account for the observed radio halo. For the 
choice of o-p = 2.3 and Xb = 0.01 the maximum of the CRp-to- 
thermal energy density XcRp is smaller than 3% for the entire range 
of the radio halo. Conservative choices for CRp spectral breaks 
have been adopted by means of Eq. <2.14t : We assume a high- 
momentu m break of ~ 2 x 10^ GeVc"' being derived from CRp 
diffusion ( Berezinskv et alJl997^ while the lower momentum cut- 
off assumes the CRp to be accelerated from a ther mal Maxwellia n 
distribution, p, ~ 3^j7m^kl\~ = 0.01 GeVc"' iMiniati""200l|). 
This choice of the lower cutoff represents the energetically tightest 
constraint because taking into account Coulomb losses would only 
weaken the energetic requirements. Moving away from the mini- 
mum energy solution, especially in the inner parts of the cluster, 

^ Though, some part of this apparent increase is an ailifact owing to 
azimuthally averaging the non-centrally symmetric synchrotron bright- 
ness distribution (for a more detailed discussion on this topic, see 
IPfrommer & EnBlij2004l) . 
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the presented energetic considerations show that the hadronic sce- 
nario is a viable explanation of the Coma radio halo as long as the 
spatially constant CRp spectral index is between 1.4 < o-p < 2.8. 

Giovannini et al. ( 1993) found a strong radial spectral steepen- 
ing from ay = 0.8 - 1.8 which would translate within the hadronic 
scenario into a CRp spectral index steepening of o-p = 1.6 - 3.6. 
If the strong steepening will be confirmed, the hadronic scenario 
will face serious challenges even when including conservative CRp 
spectral breaks. However, an absent or weaker steepening of the 
CRp spectral index e.g. from = 2.3 in the cluster centre to 
ttp = 2.8 at the outskirts of the radio halo would only double 
the CRp energy density require d to explain the r adio halo in the 
hadronic scenario. The studies o f lGiovannini et alj ri993) are based 
on two synthesis aperture radio maps obtained with different ra- 
dio telescopes. The technique of interferometric radio observations 
generally suffers from missing short-baseline information leading 
to an uncertainty of emission from larger structures: the so-called 
"missing zero spacing"-problem. This uncertainty of the surface 
brightness distribution at a single frequency is even increased for 
spatial distributions of the spectral index which represent a ratio of 
surface brightness distributions yielding to possible observational 
artefacts at the outskirts of the radio halo. Thus, future observations 
are required to decide whether the strong spectral steepening as a 
function of radius is an observational artifact or a real characteristic 
of the radio halo. 

Figure0compares radio synchrotron profiles of the Coma ra- 
dio hal£_b^;_^Gwoni_etdj_ tZfifii) which is based on observations 
bv lVenturietalJi 19901) using a sy nthesis aperture telescope with 
observations bv IPeiss et all ^1997^ using a s ingle-dish telescope. 
The statistical variance given bv lGovoni et alj 12001*) represents the 
rms scatter within concentric annuli (shown in light grey) which is 
composed of measurement uncertainties and non- sphericity of the 
underlying radio profile. Rescaling with the square root of the num- 
ber of independent beams within concentric annuli yields statistical 
uncertainties (black) however without taking into account system- 
ati cs.^ The top panel shows an comparison of th e different profiles 
bv lGovoni et al. 1.2001.) and Deiss et alj ^93) (rescaled from its 
original observational frequency v = 1.4 GHz to 326 MHz using 
the synchrotron spectral index or,, = 1.15). There seems to be an 
indica tion of spectral steepening as reported by Giovannini et al. 
il993l) . However, for simplicity we use a spatially constant CRp 
spectral index. Because of the more extended profile of the radio 
halo of the sin gle-dish observatio n, we decided to adopt the pro- 
file obtained bv lPeiss et aljfl997l) in our analysis shown in Fig.|2| 
Nevertheless, the lower panel of Fig.|4|shows a comparison of the 
energetically favoured magnetic-to-thermal energy density Xg^^^^ir) 
as a function of deprojected radius within the hadronic minimum 
energy criterion for both data sets. The tolerance regions of Xg^,-^ 
are drawn light shaded while the overlap of Xg^^.^ using the dif- 
ferent synchrotron profiles is shown dark shaded. There is no sig- 
nificant difference within the allowed tolerance regions. Together 
with the previous considerations about spectral steepening, this in- 
dicates that a moderate radially dependent spectral index does not 
significantly modify the hadronic minimum energy condition while 
a strong steepening would challenge the hadronic scenario. 



^ Due to the interferometric nature of the measurement and due to the non- 
completely synthesised aperture, we expect the true eiTor bars to be larger 
than the estimates given here suggest. However, a detailed discussion of this 
topic is beyond the scope of this work. 
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Figure 4. Top panel: Azimuthally averaged radio brightness profile of the 
radio halo in the Coma cluster as a function of imp act parameter rj_. Shown 
are the radio data at 326 MHz iOovoni et all200lh in combination with the 
Icr-error bars (black) and the surface brightness fluctu ations within concen - 
tric annuli (grey) which are based on observations by'Venturi et al.' '199(f). 
Also presented is the model profile acc ording tOtO ovoni et al. ( 2001) (dot- 
ted) and the model profile according to lPeiss etalj 119971) (dashed) which 
is rescaled to 326 MHz using the synchrotron spectral index o-y = 1.15 (c.f. 
Tablein. Lower panel: Profiles of the magnetic-to-thermal energy density 
^^rnin ('^) "'^ Coma cluster as a function of deprojected radius are shown 
within the hadronic mi nimum energy crit erion. A compaiison of the used 
synchrotron Drofiles_bv|Deisset alj_}l 997ft (solid, tolerance region medium 
grey shaded) and iGovon^^^TSMi Tdotted. tolerance region light grey 
shaded) shows no significant difference within the allowed logarithmic tol- 
erance regions (overlap is shown dark shaded). 



6 MINIMUM ENERGY CRITERIA IN A NUTSHELL 

This section provides self-consistent recipes for applying the clas- 
sical and hadronic minimum energy criterion in typical observa- 
tional situations. We present formulae for inferring magnetic field 
strengths solely as a function of observed flux per frequency, T'y, 
luminosity distance to the galaxy cluster, D, extent of the cluster 
measured in core radius, r^, observed frequency, v, and spectral in- 
dex of the diffuse synchrotron emission, Oy, where the emissivity 
scales as jy oc v""''. 

The omnidirectional (i.e. integrated over 47t solid angle) lu- 
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minosity per frequency is given by the volume integral of the syn- 
chrotron emissivity, jy. 



(6.1) 



£y = 4n \ dV jy. 

We choose a reference luminosity -£r(i+;),o = 4nD- (I + z)"' Ty^o 
which corresponds to a flux at v = 1 GHz of "Fv.o = 1 Jy for a source 
at a luminosity distance of D = 100 hjg Mpc. This corresponds 
to a cluster like Coma which is characterised by a core radius of 
~ 300 /iyg kpc. 



6.1 Classical minimum energy criterion in a nutsliell 

Applying the classical minimum energy criterion, we infer an opti- 
mal magnetic field strength by rewriting Eq. <3.1H . 



1 + £y 



2 £y,o \ rc. 



UGHz/ lO.lGeV/ 



l/(a,.+3) 



(6.2) 



Here B'^^'^^^^iay) is given by Table |5| kp denotes the ratio between 
CRp and CRe energy densities, E[ denotes the lower cutoff of the 
CRe population, and fg denotes the filling factor of the magnetic 
field in the volume occupied by CRe, which is thought to be of 
order unity. While deriving Eq. <6.2> we implicitly assumed that 
E2 » El. We also applied a lower Coulom b cutoff to the CRe 
distribution of Ei =0.1 GeV (as suggested bv lSarazinligggl) . and 
adopt a conservative choice for the cosmic ray energy scaling of 

In the case of the classical minimum energy criterion the tol- 
erance region of the magnetic field strength is given by Eq. <3.16> . 
The following substitutions might be useful when computing the 
tolerance levels of the magnetic field, crins = cr|„ ,-^/2, which are 
equally spaced in logarithmic units of B. In linear representation 
of B, this definition explicitly implies the tolerance levels given 
by exp(ln6 ± cri^g). The scaled synchrotron index is given by 
6 = (ay - l)/2, while the dimensionless magnetic energy density 



and the constant Cciass(Q'v) are denoted by 



gel 



Be 



gel 



31 i^G 



\lGHz/ 



(6.3) 
(6.4) 



Celass (o-y) = /elass (^y) /jj ' ( — ) ( ■ ^„ ) , 

-Cy.O \ ?-e.O / \ 1 GHz / 

where /eiass(Q'y) is given by Table|2| 

6.2 Hadronic minimum energy criterion in a nutsliell 

In the case of the hadronic scenario the energetically favoured mag- 
netic field strength B^^^^g is given by 

1^.0 = V^Sc = 31 nG (y^)' (6.5) 

where xb^^^ is given by Eqs. <3.28> through <3.31> . is spec- 

ified in Table |2| for a few spectral indices o-y where we assumed 
no cutoff of the CRp distribution. Provided the synchrotron index 
ffy ^ 1, there must be an upper cutoff of the CRp distribution in 
order to ensure a non-divergent CRp energy density. This might be 
obtained by means of Eq. <2.14> . Owing to the non-analytic struc- 
ture of the hadronic minimum energy criterion <3.27> in -Tb„,.„, we 
were forced to carry out an asymptotic expansion for xb^^.^ which 



Table 2. Useful numerical values for particular choices of the synchrotron 
spectral index a,, in the framework of the minimum energy criterion in a 
nutshell are given for both scenarios. Note, that a priori, we assume no cut- 
off in the CRp distribution. For spectral indices tty ^ 1 within the hadronic 
scenario, an upper cutoff needs to be introduced for deducing an equiparti- 
tion magnetic field strength in order to meet the requirements of the regu- 
larity conditions. 



/elass 



[^G] 



ghadr 

min,0 

[^G] 



/had 



0.65 


2.3 


0.6 


3.3 X ID-' 






0.75 


2.5 


0.7 


3.6 X 10-' 






0.85 


2.7 


0.8 


4.7 X 10"'' 






0.95 


2.9 


1.0 


6.4 X IQ-' 






1.05 


3.1 


1.2 


9.2 X lO-*^ 


4.0 


2.1 X 10^2 


1.15 


3.3 


1.4 


1.3 X IQ-'^ 


4.2 


1.9 X 10"^ 


1.25 


3.5 


1.7 


2.0 X 10-« 


5.3 


3.4 X 10-2 


1.35 


3.7 


2.0 


3.0 X lO"*^ 


7.3 


8.2 X 10-2 


1.45 


3.9 


2.4 


4.5 X 10"'^ 


13.0 


3.9 X 10-' 



does not admit a comparable simple scaling of the magnetic field 
as in the classical case <6.2t . 

The theoretically expected tolerance levels of the magnetic 
field, a-]„B = o"1ii.ib/2, are given by Eq. J3.32t while neglecting the 
radial dependence in the nutshell approach. The following substitu- 
tions might be useful when computing the magnetic field strength 
^nSiO corresponding tolerance region. The scaled syn- 



chrotron index is given by ^ = (ay - 
energy density of the CMB, xcmb> 
denoted by 



1)/2, while the dimensionless 
and the constant Chadr(c^v) are 



Chadr(Q'y) 



f>2 

^CMB 



„^)/^-.AfzLrf (6.7) 

" AoUcoj W^ol llGHz/ 



where /hadi(Q'y) is given by Table|2| and «e,o 



7 CONCLUSIONS 

We investigated the minimum energy criterion of radio synchrotron 
emission in order to estimate the energy density of magnetic fields 
with the main focus on the underlying physical scenario. The clas- 
sical scenario might find application for cosmic ray electrons (CRe) 
originating either from primary shock acceleration or in-situ reac- 
celeration processes while the hadronic model assumes a scenario 
of inelastic cosmic ray proton (CRp) interactions with the ambient 
gas of the intra-cluster medium (ICM) and thus leads to extended 
diffuse synchrotron and y-ray emission. 

Generally, the hadronic minimum energy estimates allow test- 
ing the hadronic model for extended radio synchrotron emission 
in clusters of galaxies. If it turns out that the required minimum 
non-thermal energy densities are too large compared to the thermal 
energy density, the hadronic scenario has to face serious challenges. 
For the classical minimum energy estimate, such a comparison can 
yield constraints on the accessible parameter space spanned by the 
lower energy cutoff of the CRe population or the unknown contri- 
bution of CRp to the non-thermal energy density. 

For the first time we examine the localisation of the predicted 
minimum energy densities and provide a measure of the theoreti- 
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cally expected tolerance regions of these energetically favoured en- 
ergy densities. The tolerance regions of the particular energy den- 
sities inferred from the classical minimum energy criterion are ap- 
proximately constant for varying magnetic field strength sg. On the 
contrary, the hadronic minimum energy criterion predicts constant 
energy densities for varying magnetic field strength in the case of 
low Eg compared to Scmb^ vvhile the tolerance region of the CRp 
energy density decreases at the same rate as the tolerance region of 
Sg increases for high Eg. 

Future observations should shed light on the hypothetical re- 
alisation of such an optimal distribution of energy densities in Na- 
ture: Combining upper limits on the inverse Compton (IC) scat- 
tering of cosmic microwave background photons off CRe within 
the ICM provides lower limits on the magnetic field strength. Un- 
ambiguous detection of the /r^-decay induced y-ray emission ow- 
ing to hadronic CRp interactions in the ICM together with the ob- 
served radio synchrotron emission yields strong upper limits on the 
magnetic field strength. These are only upper limits because the 
inevitably accompanying hadronically generated CRe could have 
a non-hadronic counterpart CRe population which also contributes 
to the observed synchrotron emission. A combination of IC detec- 
tion in hard X-rays, radio synchrotron emission, and hadronically 
induced y-ray emission therefore simultaneously enables the deter- 
mination of the CRp population as well as a bracketing of the total 
magnetic field strength and the CRe population. Applying the ap- 
propriate minium energy arguments would yield information about 
both the dynamical state as well as the fragmentation of the spatial 
distribution of the magnetic field. 

Requiring the sum of cosmic ray and the magnetic field en- 
ergy densities to be minimal for the observed synchrotron emission 
of the radio halo of the Coma cluster and the radio mini-halo of 
the Perseus cluster yields interesting results: Within the theoreti- 
cally expected tolerance regions, equipartition is possible between 
the energy densities of CRp and magnetic fields, i.e. the minimum 
energy criterion always seems to choose equipartition to be a quasi- 
optimal case. Applying the hadronic minimum energy criterion to 
the diffuse synchrotron emission of the Coma cluster yields a cen- 
tral magnetic field strength of i?coma = 2.4;^["g |J.G while in the case 
of the cool core cluster Perseus we obtain fipcseus = 8.8^5'^** \iG. 
These values agree with magnetic field strengths inferred from 
Faraday rotation which range in the case of clusters without cool 
cores within [3 \iG, 6 \iG] while cool core clusters yield values of 
~ 12 ^G iVogt & EnBlin 2003). Within the hadronic model for the 
radio mini-halo in the Perseus cluster, this results in a confinement 
for the CRp energy density of 2% ± 1% of the thermal energy den- 
sity while the magnetic energy density reaches only 0.4% of the 
thermal energy density within large uncertainties. These energetic 
considerations show that the hadronic scenario is a very attractive 
explanation of cluster radio mini-halos. 

In order to account for the radio halo of Coma in the hadronic 
scenario, the product of ScRp and Eg needs to increase by nearly two 
orders of magnitude relative to the square of the thermal energy 
density Cth towards the outskirts of the halo. Moving away from 
the minimum energy solution and adopting for instance a constant 
magnetic-to-thermal energy density, it is energetically possible to 
explain the observed synchrotron emission hadronically by only re- 
quiring the magnetic and CRp energy density to be a few per cent 
relative to the thermal energy density (and even less for the CRp in 
the cluster centre, provided Op ~ 2.3 and the cluster is isothermal). 
Such a magnetic energy density corresponds to a central magnetic 
field strength of 6 \iG. Assuming a lower magnetic field strength 
of 3 \iG corresponding to a magnetic-to-thermal energy density of 



approximately 0.5% requires the CRp energy density to be lower 
than 10% for the entire range of the radio halo. 

The considered hadronic scenario assumes a CRp spectral in- 
dex which is independent of position and thus the radio emission 
does not show any spatial variations over the clusters volume. In 
principle, one could allow for radial spectral variations of the CRp 
and thereby for the radio emission by adopting a particular his- 
tory of this population. For instance, one possible scenario would 
be given by continuous in-situ acceleration of CRp via resonant 
pitch angle scattering by turbulent Alfven waves. We discuss that 
a moderate radial steepening would not significantly modify the 
hadronic minimum energy co ndition while a confirm ation of the 
strong steepening reported bv lGiovannini et all i 19931) would seri- 
ously challenge the hadronic scenario. 

As a caveat, it should be stressed that the inferred values for 
the particular energy densities only represent the energetically least 
expensive radio synchrotron emission model possible for a given 
physically motivated scenario. This minimum is not necessarily re- 
alised in Nature. Nevertheless, our minimum energy estimates are 
also interesting in a dynamical respect: Should the hadronic sce- 
nario of extended radio synchrotron emission be confirmed, the 
minimum energy estimates allow testing for the realisation of the 
minimum energy state for a given independent measurement of the 
magnetic field strength. Within the tolerance regions, our minimum 
energy estimates in Perseus and Coma agree well with magnetic 
field strengths inferred from Faraday rotation. Under the hypothe- 
ses of correctness of the hadronic scenario, such a possible realisa- 
tion of the minimum energy state would seek an explanation of a 
first principle enforcing this extremal value to be realised in Nature. 
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